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I review the current understanding of the chiral phase transition in QCD, with particular 
emphasis on recent results in the instanton liquid model. 
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lO , 1. Introduction 

An analysis of the spectra of particles produced in heavy ion collisions at CERN and 
AGS indicates that the excited matter created in these reactions spends most of its life 
time close to the QCD phase transition Under these circumstances, we cannot expect 
that the phenomena observed in these collisions can be understood in terms of a per- 
turbative plasma of quarks and gluons occupying the reaction zone. Instead, we have to 
address the non-perturbative aspects of QCD associated with the phase transition itself. 
Furthermore, it is these phenomena that can really teach us something about the struc- 
ture of the QCD vacuum at zero temperature and baryon density. In this contribution, I 
will try to summarize our current theoretical understanding of the QCD phase transition. 
In particular, I will report on recent progress in understanding the chiral phase transition 
in the instanton liquid model. 

Lattice results on the QCD phase transition are discussed in E. Laermann's contri- 
bution, and I will not go into detail here. Nevertheless, I would like to emphasize one 
important point. QCD has a very rich phase structure as a function of the number of 
colors, the number of flavors and their masses. In particular, pure gauge QCD has a 
first order deconfinement transition, while QCD with Nf = 3 massless flavors has a first 
order chiral transition, see figure [I]. These two transitions are not connected. When the 
mass of the fermions is varied from m = to m = oo (corresponding to the pure gauge 
theory), the two first order transitions are separated by a region in the phase diagram in 
which there is no true phase transition, just a rapid crossover. Indeed, the order of the 
phase transition for real QCD, with two light and one intermediate mass flavor is still not 
established with certainty. 

The distinction between the pure gauge deconfinement and the light quark chiral phase 
transition is not just purely academic. In fact, the two transitions have completely dif- 
ferent energy scales. The chiral phase transition takes place at T c ~ 150 MeV, while 
the purge gauge transition occurs at T c ~ 260 MeV (where the scale is set by the rho 
meson mass or the string tension). In terms of energy density (and tax dollars!), this is 
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Figure 1. Schematic phase diagram for 
QCD in the Tn u ^—m s plane, from 0. The 
points show the type of transition found in 
lattice simulations with Wilson fermions. 
Note that QCD (*) appears to be in the 
first order region, while earlier simulations 
with staggered fermions placed QCD in 
the crossover region ||. 



an order-of-magnitude difference. Also, the latent heat associated with the pure gauge 
transition is rather large, on the order of 1.5GeV/fm 3 , while the chiral decondensation 
energy is Ae ~ 250 MeV/fm 3 . This has important consequences for the non-perturbative 
gluon condensate. While most of the gluon condensate is removed across the pure gauge 
transition, there is evidence that a significant part of it remains above the chiral transition. 

2. Vacuum engineering 

In Monterey T.D. Lee reminded us that the ultimate goal of relativistic heavy ion 
collisions is vacuum engineering, the removal of the quark and gluon condensates present 
in the T = \x = vacuum. Here in Germany, vacuum engineering has a long tradition. 
More than three hundred years ago Otto v. Guericke, after inventing a suitable pump, 
demonstrated the existence of air pressure by evacuating a pair of hollow semi-spheres. In 
QCD, we have to overcome the non-perturbative vacuum pressure in order to produce a 
perturbative vacuum state. The vacuum pressure is determined by the the trace anomaly 

P=~\(T^) = 7^{atsG 2 )-^Y,mf{q f q f ), (1) 

where b = HN c /3 — 2Nf/3 is the first coefficient of the beta function. Using the canonical 
values of the condensates, this relation gives p = 500MeV/fm 3 . At low temperature, the 
T-dependence of the condensates is determined by chiral perturbation theory 0] 




To leading order, the T-dependence of the quark condensate has a very simple interpreta- 
tion in terms of the number of thermal pions times the pion matrix element of the quark 
condensate, (vr|gg|7r) = —(qq)/f% 0. This means that pions act as a vacuum cleaner for 
the chiral condensate, each thermal pion removes ~ 5 qq pairs. On the other hand, the 
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gluon content of the pion is rather small. Naively extrapolating these results to larger 
temperatures, one expects chiral symmetry restoration to occur at T ~ 260 MeV, while 
the gluon condensate is essentially T-independent. This estimate is not strongly affected 
by higher loop corrections in ChPTh. Clearly, chiral perturbation theory was not meant to 
be used near the transition. Nevertheless, the result indicates that something more than 
just pions is needed to restore the symmetry at the expected temperature. It also shows 
that even if the chiral expansion is apparently convergent, the neglected (exponentially 
small) terms need not be small. 

The thermodynamics of the phase transition is usually described in terms of a bag 
model equation of state. This EOS directly incorporates the idea that the transition 
takes place as soon as the perturbative pressure from quarks and gluons can overcome the 
non-perturbative bag pressure in the QCD vacuum. For Nf = 2 and B = 500MeV//m 3 
from QTp this gives the estimate T c = [(90 j B)/(37tt 2 )] 1 / 4 ~ 180 MeV. Analyzing lattice 
thermodynamics in more detail one finds that this estimate is too large, because only 
~ 1/2 of the bag pressure is removed across the phase transition f^H. 

An even simpler approach to estimate the critical temperature is based on the idea that 
the transition occurs when thermal hadrons begin to overlap. The density of hadrons 
becomes of order 1 fm -3 near T ~ 200 MeV. This is reasonably close to T c (although, if 
T c is really 150 MeV then nh a d(T c ) ~ 0.15 fm -3 ). Nevertheless, this kind of argument is 
not correct in general. A good example is pure gauge QCD, where T c ~ 250 MeV, but 
the lightest state is at ~ 1.7 GeV, so the particle density below T c is small, on the order 
of n ~ 0.005 fm -3 . The lesson is that a first order transition occurs if the high T (QGP) 
and low T (hadronic) phases have the same free energy; the phase transition point cannot 
be inferred from looking at the low temperature phase only. 



3. QCD near T c 

Chiral perturbation theory is based on a non-linear effective lagrangian in which the 
a, the chiral partner of the pion, is eliminated. Near T c , this description is not expected 
to be useful. However, in the vicinity of a second order phase transition, universality 
implies that critical phenomena are governed by an effective Landau-Ginzburg action for 
the chiral order parameter. In the case of QCD with two (massless) flavors the order 
parameter is a four- vector 0° = (cr, 7?). Universality makes definite predictions for the 
critical behavior of (qq), the chiral susceptibility and the specific heat near T c pTTfj . At 
the moment, these predictions appear to be in agreement with lattice gauge results |I2"||, 



but the issue has not been completely settled [fL3|| . 

I would like to make a few comments concerning the role of universality arguments. 
First, it is important to clearly distinguish between the low energy chiral lagrangian (or 



the linear a- model used as an effective lagrangian, see e.g. |T3]) and the effective action 
for the order parameter near T c . The Landau-Ginzburg action is a three dimensional 
action for static modes only. It is applicable only near T c . In particular, the parameters 
in the effective action, the tt, a masses and couplings, are completely independent of 
the parameters used in the linear cr-model at T = 0. My second point concerns the 
thermodynamics of the phase transition. The effective action describes the singular part of 
the free energy only. In QCD we expect a large change in the free energy that corresponds 
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to the release of 37 (quark and gluon) degrees of freedom. This means that in practice, 
the non-universal, regular, part of the free energy will most likely dominate the universal, 
singular, contribution. 

Universality predicts the behavior of three dimensional (screening) correlation functions 
near T c . The corresponding screening masses have also been studied in some detail on the 
lattice, see section 5. In practice, however, we are more interested in dynamical (temporal) 
masses, corresponding to poles of the spectral function in energy, not momentum. These 
quantities are hard to extract on the lattice, although some exploratory studies have been 



made fI5 |. In addition to that, we have made significant progress in studying temporal 
correlation functions in the instanton liquid model (see below). The only general approach 
to the problem that we have available at the moment are QCD sum rules. 

The general strategy is easily explained. It is based on matching phenomenological 
information contained in hadronic spectral function with perturbative QCD, using the 
operator product expansion. At finite temperature, the sum rules are of the type 



c loga; 2 + c 1 (O 1 )i J + c 2 (O 2 )^ + ...= fdu 2 ^^- 2 , (4) 



where p(w 2 ) is the spectral function at p = 0, q are temperature independent coefficients 
that can be calculated in perturbative QCD and (Oj) are temperature dependent conden- 
sates. If there is a range in energies in which both the OPE has reasonable accuracy and 
the spectral representation is dominated by the ground state, we can use the sum rules 
to make predictions about ground state properties. 

In practice, this is a difficult game, even at zero temperature. At T ^ 0, additional 
problems arise because we do not know the T-dependence of the condensates and there is 
little phenomenological guidance concerning the form of the spectral functions. For this 
reason, reliable predictions can only be made at small temperature. The most systematic 



studies have been made in the vector meson channels p and a\ [|i6Hl8|. To order T 2 , 
there is no shift in the resonance masses. The only effect is mixing between the p and a\ 
channels, which is caused by scattering off thermal pions. At order T 4 , masses start to 
drop. It is interesting to note that at this order, the mass shift is not controlled by the 
quark condensate, but the energy momentum tensor in a thermal pion gas. 



4. The instanton liquid at finite temperature 

In order to make progress we need a more detailed picture of the chiral phase transition. 
In particular, we would like to understand the mechanism of the transition and the be- 
havior of the condensates and hadronic correlation functions in the transition region. In 
the following, I will argue that important progress in this direction has been made in the 
context of the instanton liquid model. In this model, chiral symmetry breaking is caused 
by the derealization of quark zero modes associated with instantons. Chiral restoration 
takes place when instantons and antiinstantons form molecules, the quark modes become 
localized and the quark condensate is zero. 

The essential assumption underlying the instanton model is that the (euclidean) QCD 
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partition function 

N f 

Z = f DA fl exp(-S)]\det(ift + im f ), (5) 

/ 

is dominated by classical gauge configurations called instantons. Instantons describe 
tunneling events between degenerate vacua. As usual, tunneling lowers the ground state 
energy. This is why instantons contribute to the vacuum energy density and pressure in 
the QCD vacuum. The instanton solution is characterized by 12 parameters, position (4), 
color orientation (7) and size (1). An ensemble of interacting instantons is described by 



the partition function 

1 - N++N- N f 

Z = amaH" / n [d A ZidUidpi d(fti)] exp(-Si n t) jjdet(i|) + im f ) , (6) 

where N + and iV_ are the numbers of instantons and antiinstantons and d(p) is the semi- 
classical instanton density calculated by 't Hooft. There are two important pieces of 
evidence that suggest that instantons play an important role in the QCD vacuum. One is 
provided by extensive calculations of hadronic correlation functions in the instanton liquid 
model jnj-|2Tfl. These correlators agree both with phenomenological information [^] and 



lattice calculations |[23|| . The second comes from direct studies of instantons on the lattice. 
An example is shown in figure Using a procedure called "cooling" one can relax any 
given gauge field configuration to the closest classical component of the QCD vacuum. 
These configurations were found to be ensembles of instantons and antiinstantons. The 
MIT group concludes that the instanton density is (N/V) ~ (1.4 — 1.6) fm~ 4 while the 
typical size is about p ~ 0.35 fm ||24|| . These numbers are in good agreement with the 
instanton liquid parameters (N/V) = lfm -4 , p = 1/3 fm proposed by Shuryak a long 
time ago p5[ . What is even more important is that hadronic correlation functions remain 
practically unchanged during the cooling process. This suggests that instantons play a 
dominant role in generating the spectrum of light hadrons. 

Recently, the role of instantons at finite temperature has been reevaluated. The semi- 
classical expression for the instanton density at finite temperature contains a suppression 
factor ~ exp(— (2N c /3 + Nf/3)(7rpT) 2 ) |2(| which comes mostly from Debye screening of 



the instanton field. From this expression we would expect that instantons are significantly 
suppressed (~ 0.2) near T c . This is not correct. The perturbative expression for the Debye 
mass is not applicable until the temperature is significantly above T c . We have already seen 
that the gluon condensate is essentially temperature independent at small T. A similar 
calculation can also be performed for the instanton density, with the same conclusion \T7\. 



This result was confirmed in a number of lattice calculations (in quenched QCD) |28|-|30| 
see figure [| The figure shows the topological susceptibility which, in quenched QCD, is 
roughly equal to the instanton density. The Pisarski-Yaffe prediction (labeled P-Y) clearly 
underpredicts the topological susceptibility near T c . The dashed curve which fits the data 
above T c corresponds to the PY- prediction with a shifted temperature T — > (T — T c ) . 

If instantons are not suppressed around T c , then the chiral phase transition has to be 
caused by a rearrangement of the instanton liquid. A mechanism for such a rearrangement, 



the formation of polarized instanton-antiinstanton molecules, was proposed in [31,32]. In 
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Figure 2. Instanton content of a typical T = gauge configuration, from ||24|| . Figs, (a) 
and (c) show the field strength while (b) and (d) show the topological charge density. The 
upper panel shows the original configuration and the lower panel the same configuration 
after 25 cooling sweeps. 



the presence of light fermions, instantons interact via the exchange of 2Nf quarks (see 
fig. ^|). The amount of correlations in the instanton ensemble depends on the competition 
between maximum entropy, which favors randomness, and minimum action, which favors 
the formation of instanton anti-instanton pairs. At low temperatures the instanton system 
is random and chiral symmetry is broken. At high temperatures the interaction in the 
spacelike direction becomes screened, whereas the periodicity of the fields in the timelike 
interaction causes the interaction in that direction to be enhanced. Schematically, the 
fermion determinant for one pair looks like 



det (ip) ~ | sin(7rTr)/ cosh(7rTr)| 



2N f 



(7) 



where r and r are the separations in the temporal and spatial direction. This interaction 
is maximal for r = and r = 0/2 — (1/2T) which is the most symmetric configuration 
of the II pair on the Matsubara torus. 

In numerical simulations, we find the transition to a correlated system in which chi- 
ral symmetry is restored at T c ~ 130 MeV |[33|| . Typical instanton configurations below 
and above T c are shown in figure |5]. The plots are projections of a four dimensional box 
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Figure 4. Schematic picture of the phase Figure 5. Instanton configurations below 
transition in the interacting instanton liquid, and above the phase transition. 



into the xr-plane. At high temperature, the imaginary time axis is short. The location of 
instantons and antiinstantons is denoted by ± signs, while the lines connecting them indi- 
cate the strength of the fermionic "hopping" matrix elements. Below T c , there is no clear 
pattern. Instantons are either isolated or part of larger clusters. Following the hopping 
matrix elements, quarks can propagate over large distances and form a condensate. Above 
T c , instantons are bound into pairs. The propagation of quarks in the spatial direction is 
suppressed and no condensate is formed. 

More details are provided in figure [| The quark condensate is practically T-independent 
at small temperature, but drops rapidly above T = 100 MeV. At small T, the instan- 
ton density rises slightlyQ It drops above T = 100 MeV, but retains about half of its 
T = value near T c . Similarly, the instanton related free energy^ does not vanish at T c . 
Instantons still contribute to the energy density and pressure above T c . 

1 This might very well be an artefact. The T = point is not a true zero temperature calculation. 

2 Figure only shows the instanton-related part of the free energy, the full free energy F = —p has to 

be a monotonically decreasing function of T. 
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Figure 6. Instanton density, free energy and Figure 7. Screening masses in the instan- 
quark condensate as a function of T. ton liquid model. 



It is interesting to study the phase diagram in more detail. For realistic QCD, the 
transition appears to be weakly first order (the discontinuity in the free energy is smoothed 
out in a finite volume and cannot be seen in figure ||). In the case of two flavors, we see a 
second order phase transition with critical exponents consistent with the 0(4) universality 
class. For more flavors, the transition temperature drops until (around Nf = 5 massless 
flavors) chiral symmetry is restored in the ground state even at T = 0. 

5. Hot hadrons 

We have studied both temporal (related to the spectral function in energy) and spatial 
(related to screening masses) correlation functions at finite temperature pi] . Figure ^ 
shows the spectrum of spacelike screening masses. First, we clearly observe that chiral 
symmetry is restored at T ~ 130 MeV. Chiral partners, like the (7r, a) and (p, a±) become 
degenerate at T c . Second, even above T c , the scalars 7r and a are significantly lighter than 
the vectors p and a\. This result is in agreement with lattice calculations |34] , |3~5|| . The 
spectrum of screening masses is also in qualitative agreement with the predictions from 
dimensional reduction (DR) |36|- |3~8|1 , but DR fails to account for the attraction seen in 



the scalar channels. 

More interesting is the behavior of temporal correlation functions, shown in figure [8[ 
The correlation functions are normalized to free quark propagation at the same tempera- 
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ture, so all correlators start at 1 for r = 0. Below T c , the data points are denoted by open 
squares (T = 0), pentagons, etc., while the high temperature points are denoted by solid 
squares and pentagons. Above T c , the period of the correlation functions becomes short, 
P = T^ 1 ~ 1.5 fm. This means that all the spectral information is contained in a very 
short interval 0<r</3/2~ 0.7 fm. This fact makes it very difficult to extract thermal 
masses from the correlation functions (or to decide whether there are narrow poles in the 
spectral function at all). 

Again, we observe that above T c the correlation functions of chiral partners become 
equal. What is more remarkable is the fact that in the regime r < (/3/2), which is not 
(directly) affected by the periodic boundary conditions, the pion| correlation function is 
almost as large as it is at T = 0. This suggest that there is still a (n, <r)-mode above T c 
39 , 40 1 . No such effect is seen in the vector channels. The small resonance contribution 
seen in the p channel at small temperature quickly melts and above T c the correlation 
function is consistent with the propagation of two independent quarks with a small residual 
chiral mass (a mass term in the vector part of the propagator that does not violate chiral 
symmetry). 

What effect causes the resonance behavior seen in the scalar channels? At zero tem- 
perature, the instanton induced interaction between quarks is conveniently discussed in 
terms of an effective four fermion interaction^ [|4l| , f42| 



d(p)dp ■ 



(2rrp)< 



2N r 



2N r 



1 



4N, 



■(far<vvo 2 , (8) 



where d(p) denotes the density of instantons. Here, if) is an isodoublet of quark fields 
and the four vector r~ has components (r, i) with f equal to the Pauli matrices acting 
in isospace. The interaction (§) successfully explains many properties of the (T=0) QCD 
correlation functions, most importantly the strong attraction seen in the pion channel. 

The effective lagrangian (H) comes from the 2Nf zero modes associated with an in- 
dividual instanton. It is derived under the assumptions that instantons are sufficiently 
dilute and completely uncorrelated. Above T c , the collective coordinates of instantons and 
antiinstantons are no longer random, but become strongly correlated. The four fermion 
interaction induced by polarized instanton-antiinstanton molecules is given by [[3^ 



mol sym 



2NI 



(^r a 7 ^) 2 + (4>T a W ^) 2 } + —(^^f + C 8 , (9) 



where the coupling constant G is determined by the number of correlated pairs (and 
their overlap matrix element) and £g is the color octet part of the interaction. G is not 
strong enough in order to cause quarks to condense. Nevertheless, molecules produce a 
significant attractive interaction in the n and a channels. 



3 The a correlation is even larger than the 7r below T c because it receives a disconnected contribution 
from the quark condensate. 

4 Strictly speaking, a flavor antisymmetric 2Nf fermion interaction. However, for Nf = 3 and broken chiral 
symmetry (either spontaneous or explicit), we can absorb two zero modes into the instanton density. 
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A problem that has received a lot of attention recently is the fate of the U(1)a anomaly 
at finite temperature [|3-|8|. Given the large rj — ir splitting any tendency towards U(1)a 
could lead to rather dramatic effects in heavy ion collisions. Two experimental signatures 
that have been discussed are the if] /it ratio ffl] measured by the WA80 collaboration [19 



and the possibility that the rf Dalitz decay jPJ contributes to the enhancement in low 



mass dileptons seen by CERES [p0| . 

The anomaly is related to the presence of zero modes in the spectrum of the Dirac 
operator. Instantons can absorb Nf left handed quarks of different flavors and turn 
them into right handed quarks, violating axial charge by 2Nf units. This is the process 
described by the 't Hooft vertex (H). Inserting the 't Hooft interaction into the rj' corre- 
lation function splits the rj' from the pionPI. Above T c isolated instantons disappear and 



5 In QCD, flavor symmetry is broken and part of the splitting comes from the strange quark mass. 
However, in the absence of the anomaly we would expect the T) — t] 1 mixing to be almost ideal (similar to 
the p — <t> system). In this case, there is a non-strange rj which is almost degenerate with the w. 
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Figure 9. Flavor mixing in the 77 — 77' system below and above the chiral phase transition. 



instanton-antiinstanton molecules do not violate U(1)a- However, the 't Hooft operator 
can induce a tunneling event (instanton) all by itself. This can be seen as follows. If we 
keep a small current quark mass, the density of isolated instantons above T c is proportional 
to m Nf . When we calculate a U(1)a violating observable, there are Nf propagators^, each 
of which has a zero mode contributing a factor l/rrif. As a result, U(1)a is broken in the 
chiral limit rrif — > 0. 

At temperatures significantly above T c , we expect screening to reduce the strength of 
the U(1)a violating interaction. Near T c , screening is not important but chiral symmetry 
restoration affects the structure of flavor mixing in the 77 — 77' system [51|, see figure |||. 
Below T c , there is strong flavor mixing between u, d and s quarks. Flavor symmetry 
is broken, (uu) 7^ (ss), but the 77' state is almost pure singlet. Above T c , there is no 
flavor mixing between non-strange and strange quarks. As a result, the eigenstates in the 
77 — 7/ system are the non-strange and strange eta components rj^s and rjs- The anomaly 
acts only on the non-strange r]NS, so the strange rjs can become light. This effect is of 
phenomenological interest, because it might enhance strangeness production in heavy ion 
collisions. 



6. Summary 

There is substantial evidence that non-perturbative effects are important in QCD, even 
above the phase transition. This evidence comes from an analysis of lattice results for the 
equation of state, the spectrum of screening masses and temporal correlation functions 
above T c . This suggest that in order to understand the transition region we need a more 
detailed picture of the transition itself. 

We have shown that the chiral phase transition can be understood as a transition from 
a disordered instanton liquid to a correlated phase of instanton-antiinstanton molecules. 
This picture is consistent with both the observation that not all of the gluon condensate is 
removed across the phase transition and with the observed spectrum of screening masses. 
It provides interesting predictions for the behavior of hadronic modes near T c . For exam- 
ple, we suggest that (ti, <r)-like modes survive the phase transition and that the structure 

6 In the (academic) case of three massless flavors, Ua(^) violation does not affect the 77' correlation function 
above T c because the third quark in the 't Hooft vertex cannot be absorbed. In real QCD, the strange 
quark can be absorbed by a mass insertion. 
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of flavor mixing in the r\ — rf sector changes near T c . 
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